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We develop, from first principles, a general and compact formalism for predicting the electromag¬ 
netic response of a metamaterial with non-magnetic inclusions in the long wavelength limit, including 
spatial dispersion up to the second order. Specifically, by resorting to a suitable multiscale technique, 
we show that medium effective permittivity tensor and the first and second order tensors describing 
spatial dispersion can be evaluated by averaging suitable spatially rapidly-varying fields each satysif- 
ing electrostatic-like equations within the metamaterial unit cell. For metamaterials with negligible 
second-order spatial dispersion, we exploit the equivalence of first-order spatial dispersion and re¬ 
ciprocal bianisotropic electromagnetic response to deduce a simple expression for the metamaterial 
chirality tensor. Such an expression allows us to systematically analyze the effect of the composite 
spatial symmetry properties on electromagnetic chirality. We find that even if a metamaterial is 
geometrically achiral, i.e. it is indistinguishable from its mirror image, it shows pseudo-chiral-omega 
electromagnetic chirality if the rotation needed to restore the dielectric profile after the reflection is 
either a 0° or 90° rotation around an axis orthogonal to the reflection plane. These two symmetric 
situations encompass two-dimensional and one-dimensional metamaterials with chiral response. As 
an example admitting full analytical description, we discuss one-dimensional metamaterials whose 
single chirality parameter is shown to be directly related to the metamaterial dielectric profile by 
quadratures. 
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I. INTRODUCTION 

Designing the electromagnetic response of an artificial 
medium is one of the main target of modern photonics 
and metamaterial science is probably the most impor¬ 
tant research field based on such skill. Basically the de¬ 
sign is made possible by the physical fact that the elec¬ 
tromagnetic field, when traveling within a nonhomoge- 
neous medium with sub-wavelength features, is not able 
to follow its spatial rapidly-varying details so that the 
field only experiences the effect of an averaged or effec¬ 
tive medium. A number of different homogenization ap¬ 
proaches have been developed for predicting the effec¬ 
tive medium electromagnetic response and they exploit 
different and suitable approximation schemes. The sim¬ 
plest homogenization technique deals with the retrieval 
of the effective parameters from the scattering properties 
of the medium M and it is based on postulating the 
equivalence between a complex metamaterial array and 
a uniform slab of same thickness with unknown consti¬ 
tutive parameters. Another homogenization technique is 
the field averaging method which is based on the averag¬ 
ing of the electromagnetic field in a metamaterial unit cell 
[9H13] and, in analogy with the retrieval technique, it is a 
numerical method for the determination of effective pa¬ 
rameters. In addition to numerical methods, mean-field 
homegenization theories are available where the effective 
parameters are evaluated from the distribution of the 
underlying metamaterial inclusions. Examples of such 
techniques are those exploiting Lorentz [14|, Clausius- 
Mossotti jTBJ, Maxwell-Garnett (ltHl approximations, or 
based on multipolar expansion [17] and source-driven ap¬ 


proach [18] . Spatial periodicity and rapidly-varying spa¬ 
tial scales are the two basic ingredients of each meta¬ 
material homogenization approach. Two different ho¬ 
mogenization techniques are obtained by assuming one 
of these two ingredients and subsequently incorporating 
the other. Therefore, starting from the photonic-crystal 
description of the structure where the spatial periodic¬ 
ity is fully taken into account, the effective medium re¬ 
sponse can be extracted in the long wavelength regime 
I19h25| . Conversely, the spatial rapidly-varying metama¬ 
terial features allows an asymptotic multi-scale analysis 
of the sample electromagnetic response which combined 
with the array periodicity yields the effective medium re¬ 
sponse f2til - l32j . 

Even though the metamaterial inclusions pattering has 
a spatial scale much smaller than the radiation wave¬ 
length, such two scales are generally not so different 
to allow a description of the effective medium response 
only comprising effective dielectric permittivity and mag¬ 
netic permeability tensors. For this reason, the effec¬ 
tive medium generally shows an additional nonlocal re- 
spose [33l - l39j yielding spatial dispersion. It is well-known 
that the nonlocal first order contribution (i.e. contain¬ 
ing first order spatial derivatives of the electric field) 
is equivalent to a reciprocal bianisotropic response [40 ] 
whereas the second order contribution can be partially 
interpreted as a correction to the effective magnetic per¬ 
meability tensor pfj El} , a phenomenon which is known 
as artificial or optical magnetism [§, . If the effec¬ 

tive nonlocality is weak, the effective medium response 
is adequately described by reciprocal bianisotropic con¬ 
stitutive relations where the chirality tensor accounts for 
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the strength of magnetic and electric polarizations cou¬ 
pling. An efficient way for observing the effect of elec¬ 
tromagnetic chirality is considering metamaterials whose 
underlying constituents’ pattering exhibits chiral asym¬ 
metry, i.e. its mirror image cannot be superposed onto 
it, and theoretical and experimental investi gati ons have 
been performed both in three-dimensional [46h49{ and 
in two-dimensional configurations (50l - f54| . Chiral meta¬ 
materials have attracted a good deal of attention since 
they can yield giant optical activity, asymmetric trans¬ 
missio n l55l. M, repulsive Casimir force [57L backward 
waves mnur and negative refractive index [60l - l63l ]. 

In this paper we theoretically investigate the electro¬ 
magnetic response in the long wavelength regime of a 
periodic composite medium, whose inclusions are non¬ 
magnetic, and we carry out the analysis by including 
spatial dispersion up to the second order. Specifically, 
by using the ratio 77 between the composite periodicity 
and the wavelength as an asymptotic expansion param¬ 
eter, we exploit a general multiscale technique for sepa¬ 
rating the fast matter scale from the slow radiation one 
and to extract the average effect of the composite on the 
electromagnetic field. We fully develop the analysis up 
to the second order in 77 and we find that the obtained 
effective medium response has a 77° order reproducing 
the know results of standard homogenization approaches 
and 77 1 and 77 2 orders accounting for first and second or¬ 
der effective medium nonlocality, respectively. The ef¬ 
fective dielectric permittivity tensor (of rank two) aris¬ 
ing from the zeroth order is obtained by averaging over 
the metamataterial unit cell specific fast-varying fields 
satisfying electrostatic-like equations. The main result 
of the present paper is that this very simple procedure 
still works for evaluating the two tensors (of ranks three 
and four, respectively) describing effective medium non¬ 
locality up to the second order and the corresponding 
electrostatic-like equations for the suitable fast-varying 
fields to be averaged are here derived and discussed. As 
a consequence we obtain a simple scheme for evaluat¬ 
ing the effective medium response up to the second order 
which is based on solving a small number of Poisson equa¬ 
tions on the metamaterial unit cell. In particular we ob¬ 
tain a very compact expression for the three-rank tensor 
describing first order spatial nonlocality which remark¬ 
ably reveals that it can be evaluated by using the same 
fast-varying fields appearing in the zeroth order descrip¬ 
tion, thus avoiding the requirement of solving additional 
electrostatic-like equations. Therefore the ensuing pro¬ 
cedure for deducing the medium response up to the first 
order is remarkably simple and compact. For metamate¬ 
rials with negligible second-order spatial dispersion, we 
combine the obtained description with the equivalence of 
first-order spatial dispersion and reciprocal bianisotropic 
electromagnetic response to obtain a compact and sim¬ 
ple expression for the effective medium chirality tensor. 
Using such expression we investigate the impact of the 
composite patterning symmetry on electromagnetic chi¬ 
rality. Specifically we obtain that if the composite does 


not show chiral asymmetry, i.e. its mirror image can be 
superposed on it, the effective chiral tensor vanishes with 
two remarkable exceptions corresponding to the specific 
cases where the rotation needed to superpose the mir¬ 
ror image onto the structure is either a 0° or 90° rota¬ 
tion around an axis orthogonal to the reflection plane. 
In these two situations the composite is geometrically 
achiral and it nonetheless shows electromagnetic chiral¬ 
ity whose chiral tensor turns out to be that of pseudo- 
chiral-omega medium. Such general result encompasses 
in particular the relevant cases of two-dimensional and 
recently considered one-dimensional metamaterials show¬ 
ing chiral response. As an example admitting full ana¬ 
lytical description, we discuss one-dimensional metama¬ 
terials whose single chirality parameter is shown to be 
directly related to the metamaterial dielectric profile by 
quadratures. 

The paper is organized as follows. In Sec. II we dis¬ 
cuss our approach to the effective medium description of 
a periodic composite in the long wavelength regime in¬ 
cluding spatial nonlocality. This is done examining the 
general multiscale technique in Sec. II.A and applying 
it to the zeroth, first and second orders in Sections II.B, 
II.C and II.D, respectively. In Sec. Ill we sum up the 
results obtained in Sec. II and we stress the simplicity 
and compactness of the proposed approach. In Sec. IV 
we consider media whose second order contribution to 
spatial nonlocal response can be neglected and we focus 
on electromagnetic chirality. Specifically in Sec. IV.A 
we adopt the bianisotropic description of the medium 
to deduce a simple expression for the chirality tensor, 
in Sec. IV.B we examine the relation between geomet¬ 
ric and electromagnetic chirality of metamaterials and in 
Sec. IV.C we specialize such general analysis to the cases 
of two-dimensional and one-dimensional chiral metama¬ 
terials. In Sec. V we focus on one-dimensional chiral 
metamaterials and we obtain a closed-form expression 
for the single parameter ruling electromagnetic chirality 
of such media. In Sec. VI we draw our conclusions. 

II. HOMOGENIZATION THEORY 
A. Electromagnetic field multiscale analysis 

Let us consider propagation of a monochromatic elec¬ 
tromagnetic field through an unbounded metamaterial 
whose underlying non-magnetic inclusions (both metal 
and dielectric) are arranged on a lattice whose period is 
much smaller than the radiation wavelength A. The elec¬ 
tric E and magnetic H field amplitudes satisfy Maxwell 
equations 

V x E = zw/ioH, 

V x H = —icueosyE (1) 

where time dependence e~ lult has been assumed and 
the medium relative dielectric constant e r (r) is a peri¬ 
odic complex function having the same periodicity of the 
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metamaterial. Due to the rapidly varying dielectric spa¬ 
tial modulation it is convenient to introduce the param¬ 
eter r) = d/X where d is the largest of the lattice basis 
vector lengths. Following the standard multiscale tech¬ 
nique, it is possible to exploit the condition 77 -C 1 to 
develop an asymptotic analysis of electromagnetic prop¬ 
agation. Accordingly, we introduce the fast spatial coor¬ 
dinates R = r/ 77 , regarded as independent from the slow 
coordinates r. Any field /( r,R) can be decomposed as 
/(r,R) =7(r) + /(r,R) where the overline denote the 
spatial average over the metamaterial unit cell, i.e 

7(r) = i^d 3 R/(r,R) (2) 


field orders only through the averaged term eE„. Mul¬ 
tiplying each of Eqs. © for 77 " and summing over n we 
obtain 

V x E = iw/xoH, 

V x H = -iuD (7) 

where E = E «+ and H = H »*i" is the 

overall averaged electromagnetic field and 

OO 

D = £o£E + 'y ' AD n (8) 

n —0 

where 


where C is the unit cell and V is its volume scaled by 
77 s , and the tilde denote the rapidly varying zero mean 
residual, i.e. / = / — /. In our scheme the relative 
dielectric permittivity depends on the fast coordinates 
only so that we set 


e(R) = e r (? 7 R) (3) 

and it can be decomposed as e(R) = e + e(R). Follow¬ 
ing the general multiscale technique Q, we set for the 
electromagnetic field amplitudes A = E, H 


A = E 

n—0 


[A n (r) + A„(r,R) 


(4) 


or, in other words, the field amplitudes are expanded 
in powers of 77 and, for each order, the average on the 
unit cell (terms with an overbar and depending on the 
slow coordinates only) is separated from the zero mean 
residual (terms with an overtilde and depending on both 
slow and fast coordinates). After substituting Eqs.® 
into Eqs.® and noting that V —> V + 4Vr it is possible 
to extract equations for each order and, for each of them, 
to separately balance the averaged contributions and zero 
mean residuals. As a result, for the averaged equations 
we obtain 


V x E ra — iu) oH n , 

V x H ra = iuj£ q (eE n -|- eE n ) (5) 

whereas for the zero mean residual equations we have 


Vj^ x E 0 
Vr x Ho 
Vr x E„ + i 
Vr x H n+ i 


0, 

0, 


—V x E„ + iui/j, oH„ 
-V x H„ 


—iioe 0 


(e — e)E n + eE ra — eE n (6) 


AD„ = e 0 £E n 77 n . (9) 


Therefore the slowly varying electric and magnetic field 
amplitudes satisfy the macroscopic Maxwell equations 
® with the slowly varying displacement vector D of 
Eq.® which has two contributions, the former due to 
the spatial average of the dielectric profile and latter due 
to the dielectric modulation. The latter contribution is 
obtained by summing the spatial average of the rapidly 
varying fields E„ multiplied by the dielectric permittiv¬ 
ity arising from each order. Therefore, in order to obtain 
an effective medium description of the metamaterial re¬ 
sponse, such rapidly varying fields have to be related to 
the slowly varying ones through Eqs. d®). 

From the structure of Eqs. © we note that the fields 
E„ + i and H„ + i have to be evaluated recursively from 
the knowledge of E„ and H n . In order to accomplish 
this task it is convenient to deduce equations involving 
the fields divergence. After applying the operator Vr- 
to both the third and fourth of Eqs.® we obtain 


V • (V R x E n ^ 

V ■ x H„^ 


-7W/7 0 V R • H„, 


*W£qV R 


(e ~ s)E n + £E n (10) 


where we have used the identity V R ■ (V x A) = —V ■ 
(Vr x A). Setting n = 0 and using the first and the 
second of Eqs.®, Eas. ([T!Jl) become 


V r -H 0 = 0, 

Vr-(£E 0 ) = ~ (Vre) -E 0 . (11) 


Relabelling n —tn+\ in Eqs. ca and substituting the 
expressions for V R x E„+i and V R x H„ + i from the 
third and fourth of Eqs. ® we obtain 


Vr H 


n+1 


7 R ■ ( £ ®«+i) 


= -V 

= -V' 


R 4T„, 


(s ~ £)E n + £E n — eE„ 




E 


n+1 


( 12 ) 


where n = 0 , 1 , 2 , Equations ® do not contain the for n = 0,1, Equations ® together with Eas. lfTTl) 

fast coordinates R and are coupled to the rapidly varying and m can be used to evaluate the rapidly varying 
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fields of order n + 1 once those of order n are known 
and these fields are linearly dependent on the slowly- 
varying fields. Therefore, the discussed homogenization 
technique allows to obtain an effective medium descrip¬ 
tion of the metamaterial respose up to the desired order 
in ?y. In this paper we will provide an effective medium 
description up to the second order. 


B. Zeroth order 

From the second of Eqs. © and the first of Eqs. ED 
we note that the field Ho has no sources and therefore it 
vanishes. On the other hand the first of Eqs.© and the 
second of Eqs. ED state that the field Eo is conservative 
and it is produced by the spatial modulation of e and the 
slowly varying field Eo- Therefore we have 

H 0 = 0, 

E 0 = ei(difj)E 0j (13) 

where the sum is hereafter understood over repeated in¬ 
dices, e; is the unit vector along the i-th direction, di is 
the partial derivative along Xi = e,; • R, Eoj — Gj • E 0 
and the functions fj satisfy the equations 

Vr ' (^R fj) = -dje. (14) 

where j = 1, 2, 3. Note that the source term of Eci. (fl4l) 
has the same metamaterial periodicity and its spatial av¬ 
erage vanishes, i.e. dj£ = 0 , so that the fields fj are 
periodic functions with the same metamaterial period¬ 
icity and they are defined up to an arbitrary constant 
(see Appendix A). Inserting the second of Eqs. ED into 
Eq. © with n = 0we obtain 

XD 0 =s 0 e i Ae\ e / f) Eo j . (15) 

where As^j^ 1 = edifj so that the zeroth order contribu¬ 
tion to the displacement vector of the effective medium 
amounts to a tensor correction to the average dielectric 
permittity. This result agrees with Ref. (28j where the au¬ 
thors consider the homogenization of a two-dimensional 
structure up to the zeroth order in 77 . Note that, as shown 
in appendix B, the rank-two tensor Ae'/j-j satisfies the 

relation Ae = Ae and this is the correct sym¬ 
metry property of the correction to the dielectric permit¬ 
tivity tensor, in agreement with the Onsager symmetry 
principle [65| . After symmetrizing the pair of indices ij , 

i.e. by setting Ae^^ = \ (a e[j ff ' 1 + Ae ^/^, we get 

A ^tj^ = ~5ij£ + -(Qij + Qji ). (16) 

where we have set 

(17) 


for later convenience. It is worth noting that the expres¬ 
sion for only involves spatial derivatives of the fields 
fj so that, at the zeroth order, the undefined additive 
constant in fj (see Appendix A) do not affect the effec¬ 
tive correction to the permittivity tensor of Eg. (1161) . In 
addition, as a consequence of Eq. (fl4l) we have the diver¬ 
genceless condition 

diQij = dje + V R • (eV R /,) = 0. (18) 

C. First order 

In order to obtain the equations for the fast varying 
fields of the first order we use the third and fourth of 


Eqs.© and Eqs. ED with n 

= 0 thus getting 


Vr H, = 

0 , 



Vr x Hi = 

(Qij 

Qij) Eq j\ , 


V R • (eEr) = 

{Qij Qij ) 

dE{)j /rv x “T7 

(die)E u , 

OXi 

V R x 

- 

0 

(19) 

where 

we have used Eqs. El 

) and (fT71) and 

we have ex- 


~ / QE 

ploited the evident relation V x Eq = — V R x I e; fj d f 3 


Note that we use the symbol d/da ij to label the partial 
derivative operators with respect the slowly varying co¬ 
ordinates Xi in order to avoid confusion with di = d/dXi 
which we use for partial derivative operators with respect 
the fast varying coordinates. 

The first and the second of Eqs. ED provide the field 
Hi since 

Vr • (■eiQijEoj ) = (diQij) Eoj = 0, (20) 

where Ea. (fl8l) has been used, is the compatibility condi¬ 
tion of the system. Therefore we set 

= • [Vr X E 0j (21) 

XUJ /JjQ 

where the functions Ajj satisfy the magnetostatic equa¬ 
tions 

diAij — 0, 

V R Ay = —fcg (Qij — Qij ) (22) 

where ko = w/c, which have to be solved with the pres- 
critions that Ajj has the metamaterial periodicity and it 
is defined up to an arbitrary constant (see Appendix A). 
Note that if Aij satisfies the second of Eas. (l22l) then 

(d^j) = —kydiQij = 0 (23) 

where Ea. (fT51) has been used and this implies that 
diA^ = const, since the constant is the only harmonic 


Qij — £ ij T di fj ) , 








5 


function which has the metamaterial periodicity. In addi¬ 
tion diAij has vanishing average since it is the derivative 
of a periodic function and therefore diA^ = 0. Therefore 
the first of Eqs. (15121) is automatically satisfied and has not 
to be additionally required. 

The third and the fourth of Eqs. © provide the field 
Ei and it is easily seen that it can be expressed as 


Ei = e,; 


(difj) E\j + ( Si r fj A 8iW r j) 


dE 0j 


dx r 


(24) 


where the functions W r j satisfy the equation 

Vr • (eV R W rj ) = -dr (efj) - (Q rj ~ Q^) ■ (25) 


It is important noting that field Ei has, by its very def¬ 
inition, zero average and, from Ea. (l24l) this forces the 
undefined additive constant of f 3 to be set equal to zero, 
i.e. in such a way that 


fj = 0. (26) 

Note that the source term of Eci. (l25ll . in its right hand 
side, has vanishing spatial average so that the equation 
is structurally equivalent to Eq. m and therefore W r j 
has the metamaterial periodicity and it is defined up to 
an additive constant (see Appendix A). Inserting Ea. (l24l) 
into Eq.© with n = 1 we obtain 


ADi = eo bj 


rjAs^^Eij A a 


(eff) 

ijr 


dE 0j 

dx r 


(27) 


where ol>^ r E = r/e ( Si r fj A diW r j) so that the first order- 
contribution to the displacement vector of the effective 
medium provides a tensor correction to the average di¬ 
electric permittivity which has the same structure as its 
zeroth order counterpart (see Eq. (fl5l) ) and a contribution 
which is linear in the derivatives of the slowly varying ze¬ 
roth order electric field. This former contribution is re¬ 
sponsible for the first order nonlocal part of the effective 
medium dielectric response. 

The rank-three tensor otf/ r E can be conveniently 
rewritten as (see Appendix B) 

atC/ r f) = viQnfj - Qrjfi) (28) 

which reveals that the nonlocal contribution of the effec¬ 
tive medium response up to the first order can be pre¬ 
dicted directly from the knowledge of the functions fi 
thus avoiding to solve Eas. (l25l) . Note that Ea. (l28l) also 
shows that 


(e//) = (e//) 

ijr jir 


(29) 


which is the correct antisymmetric property, in agree¬ 
ment with the Onsager symmetry principle, for the third 
order tensor associated with first order spatial nonlocal¬ 
ity of the effective medium response [65]. 


D. Second order 


The equations for the fast varying electric field of the 
second order are obtained from the third of Eqs.© and 
the second of Eqs. m with n = 1 so that 


dEij ( d 15 \ d*E 0j 


' ( £ ® 2 ) — (Qij Qij ) (Pirj E irj) 

-( di£)E 2i , 


dXidXr 


V R x E 2 = V R x 


R 


e; f : 


dE U , T zr d2E 0i , a p 

' j -r- VVrjT; --T 

OXi OXiOXr 


(30) 


where we have used Eqs. m and d24J) , we have set 

Pirj =£ (Sirfj+ 8iW rj ) (31) 

and we have used the evident relation V x Ej = - V R x 
W T j jf-jr 1 -'] ■ Note that, since we are in- 


* {fj 


0 dxi r 0 dxidx r 

vestigating the effective medium response up to the sec¬ 
ond order if 2 , we are not considering the equations for 
the second order magnetic field H 2 since it provides con¬ 
tributions to the effective medium response only of order 


T 


It is easily seen that the field E 2 satisfying Eas. (l30l) is 

dE i j 


E 2 — e,- 


(difj) E 2 j A ( Sirfj A diWrj) 


dx r 


A (Aij A diVj) Egj A if is if Ay A diR sr j) 


d 2 E< 


0 j 


dx s dx r 


(32) 


where the functions R sr j and Vj satisfy the equations 

V R • = d s {sW r j ) {P sr j P sr j J , 

Vr • (sVrUj) = - di (i sAij ). (33) 


It is important noting that field E 2 has, by its very def¬ 
inition, zero average and, from Eg. (15121) this forces the 
undefined additive constants of A, 3 and W r] to be set 
equal to zero, i.e. in such a way that 


Aij — 0, 

W~ = 0. (34) 


Note that the source terms of Eas. (l33l) have vanishing 
spatial average so that such equations are structurally 
equivalent to Eq. m and therefore R sr j and Vj have the 
metamaterial periodicity and they are defined up to an 
arbitrary constant (see Appendix A). Inserting Eci. (15121) 
into Eq. © with n = 2 we obtain 


AD 2 = £oej 


jAdd/fe , „Aeff)dEij 
V A£y E 2j A va ijr 

d 2 E 0j 




0 j 




{eff) 


lJsr dx s dxr 


(35) 
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where 7 = e (Aij + diVj) and f[ e / s P = 

2 --------- 

\e ( Si S W r j + 6 ir W S j + diR sr j + diR rs j). Note that 

in the fourth term in Ea. (l35l) . PfpJ 1 has been sym¬ 
metrized with respect the indices sr since it is contracted 
with the symmetric mixed partial derivatives tensor. 
From Eg. (1351) we note that second order contribution to 
the displacement vector of the effective medium contains 
a correction to the effective dielectric permittivity tensor 
and a correction to the first order nonlocal response 
which which have the same structures as their zeroth 
and first order counterpart (see Eqs. <0 and (l27l) l. In 
addition it provides a novel contribution to the effective 
dielectric permittivity which is physically due to the 
effect of the rapidly varying first order magnetic field 
and a contribution which is linear in the mixed second 
order derivatives of the slowly varying zeroth order 
electric field, the latter being responsible for the second 
order nonlocal part of the effective medium dielectric 
response. 

The rank-two tensor and the rank-four tensor 

PijsP can rewritten as (see Appendix B) 


7 ( e //) 

hj 


o(eff) 


1 


— , 2 i^sAri) , 

Kg 
Q 

7 1 


— a (QriWsj Qsi^^rj Qrj^^si Qsj^^ri) 


- {Prsj + Psrj ) + fj (Prsi + Psv 


(36) 


so that the second order contribution to the effec¬ 
tive medium response can be predicted without solving 
Eas. (l33l) . It is worth noting that the relations Pp^ = 

and = PplP hold and they are the correct 

symmetric properties, in agreement with the Onsager 
symmetry principle, for contributions to the effective di¬ 
electric tensor and for the rank-four tensor associated 
with second order spatial nonlocality (65], respectively. 


III. EFFECTIVE MEDIUM RESPONSE UP TO 
THE SECOND ORDER 


The slowly varying electric and displacement fields up 
the second order are E = E^i = E 0 + Ei 77 + E 2 ? 7 2 and 
D = = D 0 + D 177 +D 277 2 so that, inserting Eqs. m 

m and (l35l) into Eq. ([ 8 ]) (whose series is truncated up to 
the second order) and adding suitable higher order terms 
for restoring the field Ei in the nonlocal contributions, 
we obtain 


where, from Eas. dTbl) . (l28l) and (l36l) . we have 

£i 

<0 


1 2 

[ e /^ = pQij + Qji) +-j^{d s A r i)(d s A r j), 




Op/ r f) = ViQrifj - Qrjfi) 


PPP = \{QriW sj + Q si W rj + QrjW s i + Q sj W ri ) 


- “~[fi (Prsj + Psrj) + fj (Prsi + Psri )]• 


(38) 


where 


Qij — E ( Sij T difj) , 

Pirj = E&rfj + diWrj). (39) 


In order to sum up the results obtained so far for clarity 
purposes, we here report Eq. ED. the second of Ea. (l22l) 
and Ea. C5l) for the fields fj, A i3 and W T1 necessary for 
evaluating the effective medium tensors, namely 

V R ' ( £V R fj) = 

VrAj = — A) {Qij ~ Qij ) > 

Vr • (eV R W rj ) = -dr (efj) - (Q rj - Q~) , (40) 

whose solutions have to be determined with periodic 
boundary conditions on the edges of the metamaterial 
unit cell and with vanishing spatial average. 

Equations (ED, (EED and SOD are the main result of 
this paper and they show that the effective dielectric 
tensor and the tensors and /3( e fP describ¬ 

ing first and second order spatial dispersion can easily 
be evaluated after solving Eas. (HUl) for a given periodic 
dielectric profile e(R). In turn these equations share a 
common electrostatic structure and therefore they can 
be faced through a number of well known theoretical 
and numerical schemes. In addition we note that, to 
the best of our knowledge, the second and the third of 
Eas. (l38l) contain the simplest expressions for the tensors 
of e ff) and p e fP since they are easily obtained by aver¬ 
aging suitable mesoscopic fields. We conclude that, the 
discussed homogenization scheme provides a particularly 
simple way for predicting and analyzing, in addition to 
the local dielectric response, the spatial dispersion prop¬ 
erties of an arbitrary tridimensional periodic dielectric 
medium in the long wavelength regime up to the second 
order. 


IV. ELECTROMAGNETIC CHIRALITY 
A. Reciprocal bianisotropic metamaterial response 


Di = £ 0 




Jeff) 


dJP 

8x r 


■A 


(e//) ® Ej 
ijrs dx,.dx f 


( 37 ) 


Let us focus here on the situation where the contri¬ 
butions of order if in Eg. (1371) can be neglected. It is 
well known that in this case the effective medium first 
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order nonlocal response is equivalent to a reciprocal bian- 
isotropic response which is referred to as a electromag¬ 
netic chirality floj . In order to evaluate the chiral ten¬ 
sor let us note that Maxwell equations for the effective 
medium are left invariant by the transformation 


very general and tridimensional situation it is convenient 
to express the chirality tensor as 0 

k = \tt{k)I + N + J (47) 


E' = E, 

B 7 = B, 

D 7 = D +V x V, 

H 7 = H-twV, (41) 

where V is an arbitrary vector field, since Eqs. 0, after 
defining B = /jqH, can be rewritten as 

V x E 7 = iwB 7 , 

V x H 7 = —iwD 7 . (42) 

This shows at the same time that the fields D and H are 
not uniquely defined (as opposed to the fields E and B) 
and that, for each possible vector field V, the constitutive 
relations 


tV _ A „ ( _(e//) cv i „.( e //) dE'j 
D - e t £ 0 l Eij h j + a ijr 

B 7 = fio (icuV + H 7 ) 


Vx V, 


(43) 


characterize the electromagnetic behavior of the effective 
medium. Choosing V = e* (jj^ K %jEj \ for the arbitrary 
vector field V, where 


Ka — kn I — — e 


iOL, 


(eff) 


-e a^ ef ^5- 

^nmq^nmq Vij 


(44) 


and Cijk is the Levi-Civita tensor, and exploiting the cru¬ 
cial antisymmetric property Oif^ = Eas. (El) 

become (see Appendix C) 

D 7 = £ 0 £ ,{eff) W - -K T W, 

c 

B 7 = -kES + ^oH 7 (45) 

c 

where e'C//) = + n T n. Equations (l45l) shows that 

the effective medium has a reciprol magneto-electric cou¬ 
pling or electromagnetic chirality described by chiral ten¬ 
sor k of Eg. (El) . Inserting the second of Eas. (1551) into 
Ea. (l44|) we obtain 


Kij = r)k 0 


jEfm + 


~imn&jq 


c mqn u ij 


Si, Sfmdqfn 


(46) 

This relation is among the main results of the present 
paper and to the best of our knowledge it is the simplest 
expression for the effective medium chiral tensor since it 
is directly evaluated using the functions fi obtained by 
solving the first of Eas. (HUl) . Since we are dealing with a 


where TV(re) = Ku is the trace of the chirality tensor, / 
is the identity tensor, iV is a symmetric trace-free tensor 
and J is an antisymmetric tensor. Such decomposition 
allows to classify the reciprocal bianisotropic response 
[40l | of the effective medium and the three main classes 
are Tt(k) ^ 0, N ^ 0, J = 0 for chiral media, Tr{n) = 
0,iV ^ 0, J = 0 for pseudochiral media and Tr(n) = 
0, IV = 0, J ^ 0 for omega media. 

For the chirality tensor of Eq. (Ell we obtain after some 
straightforward algebra 


Tt(k) r/ko ( ^£mqn 1 ^fm^qfni 


Nu = 


J / 7 - 


Vko 


2 (timn^jq 4“ tjmn^iq) T ^mqn&ij 

4 “ 77 (^imn^jq £jmn$iq) 


^fm^qfrn 

sfmdqfn- 

(48) 


B. Spatial symmetries and chirality tensor 
structure 


A structure is chiral if it is distinguishable from its 
mirror image or, in other words, if its mirror image can¬ 
not be rigidly superposed onto it. If the structure is not 
chiral it is necessary invariant under a set of geometri¬ 
cal symmetries which affect the structure of the chirality 
tensor k. In order to investigate the relation between 
geometrical and electromagnetic chirality it is necessary 
to identify the symmetries of the dielectric profile which 
either entail the vanishing or select a specific structure 
of the chiral tensor. 

Let us suppose that the metamaterial is achiral or, 
in other words, that the mirror image of the dielectric 
profile through a plane orthogonal to the unit vector s can 
be rigidly superposed onto the original profile through 
the composition of a rotation of an angle 9 around the 
unit vector r and a translation T. This implies that the 
dielectric permittivity is left invariant by the composition 
of these three transformations or 

e(R') = e(R) (49) 

where R 7 = A'e,; is the image of the point R = A,;e,; 
obtained through the spatial tranformation 

X[ =Ti + RijXj (50) 

where T) are the components of the translation vector 
and Rij = RitSkj (or in matrix form R = RS) is the 














improper rotation obtained by composing the reflection 
S and the rotation R whose matrices are 


Snm — &nm 2 5 n S m 

Rnm = i$nm ^ViTm) COS 6 tnkm'Y'k Sm 6. 

(51) 

In Appendix D we fully investigate the way the symmetry 
of Eci. (l49l) affects electromagnetic chirality and as a re¬ 
sult we obtain that such symmetry implies the vanishing 
of the chirality tensor apart from two specific situations 
where it imposes restrictions on the chirality tensor struc¬ 
ture. With reference to Fig.l, these cases are: (a) the ro¬ 
tation is trivial 0 = 0 or the rotation unit vector is orthog¬ 
onal to the reflection unit vector r-s = 0; (b) the rotation 
and the reflection are such that sin 2 (|) (r ■ s) 2 = i. 

In order to discuss such symmetric situations where the 
chirality tensor does not wholly vanish it is convenient to 
consider the orthonormal basis 


u« = 


i(2) = 


i(3) = 


U< 2 > X #, 


u 


(3) 


x r 


yj 1 — (r • fb 3 ))“ 

— sin (|) s x r + cos (|) s 
V^-sin 2 (f) (f -s) 2 


(52) 


Note that for 9 — 0 the rotation is the identity trans¬ 
formation so that the rotation vector r is in this case an 
arbitrary unit vector. 

In case (a), the dielectric profile is invariant under a 
pure reflection through a plane orthogonal to s or a suit¬ 
able rotation around an axis belonging to the reflection 
plane is in order for restoring the dielectric profile after 
the reflection. As it is shown in Appendix D, in this case 
the unit vectors of Eos. (l52l) are such that Ru^ 1 ) = fib 1 ), 
Rib 2 ) = ib 2 ) and Rib 3 ) = — ib 3 ) so that the improper ro¬ 
tation R is a pure reflection through the plane spanned 
by the unit vectors ib 1 ) and ib 2 ) and a rotation around 
ib 3 ) of an angle '0 = 0°. In other words case (a) deals 
with the situation where the medium admits a plane of 
reflection symmetry. In this case the structure of the chi¬ 
rality tensor allowed by this symmetry is (see Appendix 

D) 


K = v + k (31) u (3) u (1)t 

+ k^u^u^ t + K ( 32 )u (3 )ib 2 ) T ) (53) 

where «b 13 ), «b 31 ), «b 23 ), fib 32 ) are four independent com¬ 
plex scalars. From Eg. (1551) it is evident that Tr(n) = 0 
so that if the metamaterial admits a plane of reflection 
symmetry it is a pseudo-chiral-omega medium. 

In case (b) the medium profile is restored, after the 
reflection, by a specific rotation in such a way that 
sin 2 (!) (r • s) 2 = As it is shown in Appendix D, 
in this case the unit vectors of Eqs. (El are such that 





FIG. 1: (Color on-line) A reflection (dotted lines) through a 
plane orthogonal to the unit vector s followed by a rotation 
(circular arrows) of an angle 9 around the unit vector r is 
equivalent to a reflection through the plane spanned by the 
unit vectors fb 1 ) and u (2 ^ and a rotation around the unit 
vector u (3) of an angle 0 such that cos 0 = 1 — 2 sin 2 (|) (r ■ 
s) 2 . If the medium dielectric profile is invariant under the 
composition of the reflection and the rotation, three different 
situations are relevant, (a) The rotation is trivial 9 = 0 or 
r ■ s = 0 so that 0 = 0. (b) The rotation and the reflection 
are such that sin 2 (|) (r ■ s ) 2 = | so that 0 = ^ or 0 = 

(c) Conditions (a) and (b) are not satisfied. In case (c) the 
chirality tensor vanishes whereas in cases (a) and (b) it has a 
reduced structure. 


Rfb 1 ) = u^ 2 ), Ru) 2 ) = — ib 1 ) and Rib 3 ) = — tb 3 ) so that 
the improper rotation R is a reflection through the plane 
spanned by the unit vectors ib 1 ) and £b 2 ) and a rota¬ 
tion around the unit vector u (3 ) of an angle 0 = 90°. 
In other words case (b) deals with the situation where 
the medium is left invariant (modulus a translation) by 
a reflection through a plane and a 90° rotation around an 
axis orthogonal to such plane. In this case the structure 
of the chirality tensor allowed by this symmetry is 


k = r] 


cW 


(u^u 1 


(2 )T 


tb 2 )uW T 
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+«(“) (u«u« T - u^u^ T ) 


(54) 


where are two independent complex scalars. 

From Eci. (l54l) it is evident that Tr(n) = 0 so that if the 
metamaterial has the symmetry of case (b) it is pseudo- 
chiral-omega medium. 

We conclude that if the medium can be superposed to 
its mirror image the chirality tensor vanishes unless the 
rotation needed to restore the dielectric profile is around 
an axis orthogonal to the reflection plane and the rotation 
angle is either 0° or 90° and in both cases it is a pseudo- 
chiral-omega medium. 


C. Two-dimensional and one-dimensional 
electromagnetic chirality 

As an application of the above symmetry analysis let 
us consider two relevant situations where the medium 
is geometrically achiral and has a chiral electromagnetic 
response nonetheless. 

The first situation is that of a medium which is invari¬ 
ant under translations along an axis, say the Z— axis, so 
that its underlying dielectric profile is e = s(X,Y). This 
implies that the structure is left invariant by a reflection 
though the plane Z = 0 followed by a rotation of zero 
angle 9 = 0 around any axis orthogonal to the Z- axis, so 
that we set s = e z and r = e^. Using such unit vectors, 
the basis vectors of Eas. (l52l) becomes 

= •> 

= &y> 

= e*. (55) 

so that since the medium symmetry is of the kind (a) 
(see the last paragraph), from Eas. (l53l) we have 

0 0 rjK ( 13 ) \ 

0 0 77«( 23 ) ] (56) 

rf K ( 31 ) t y/A 32 ) 0 ) 



',( 2 ) 


U 


(3) 


so that, in agreement with the considerations of the last 
paragraph, the metamaterial behaves as a pseudo-chiral- 
omega medium. Since such geometrically achiral medium 
is translationally invariant we conclude that Eq. (1561) con¬ 
tains the most general expression for the chirality ten¬ 
sor of a two-dimensional medium, i.e. it describes two- 
dimensional electromagnetic chirality. This situation also 
encompasses the case of a planar metamaterial [50i - l54l | 
(i.e. a very thin metamaterial sheet) since, after the re¬ 
flection through a plane parallel to the medium plane, 
the sheet can be made to superpose the original sheet 
through a zero-angle rotation and a translation (which 
does not affect the chiral tensor structure) along the axis 
normal to the medium plane. 


The second situation is that of a medium which is in¬ 
variant under rotations around an axis, say the Z— axis, 
so that its underlying dielectric profile is s = e(i?, Z) 
where R = \/ X 2 + Y 2 . This implies that the structure 
is left invariant by a reflection though any plane contain¬ 
ing the Z axis followed by a rotation of zero angle 0 = 0° 
around the Z- axis. Therefore we set s = cos ^e^+sin <pe y , 
where <fr is an arbitrary angle, and r = e z . Using such 
unit vectors, the basis vectors of Eqs. (El becomes 


',(2) 


i(3) 


e*, 

sin (f>e x — cos <fre y , 
cosc/)e x + sin 0e y , 


(57) 


so that since the medium symmetry is of the kind (a) 
(see the last paragraph), from Eqs. El we have 


k = rj 


(«( 23 ) + k;( 32 ) ) sin (f> cos <j> 
—r( 23 ) cos 2 < j > + k ^ 32 - 1 sin 2 i j > 
/J 13 ) cos <j) 


K ( 23 ) gin 2 (j) _ k ( 32) cog 2 ^ K (3i) cos (jj 

— (k/ 23 ) -|- r( 32 ) ) sin <fi cos 4> K*- 31 ) sin</> 
K ( 13 ) s in</> 0 


(58) 


where K^ 13 \<f>), /c^ 31 -* hS 23 ^^) and K^ 32 \(f>) are func¬ 
tions of (j>. The entries of the chirality tensor of Ea. (l58l) 
cannot depend on the arbitrary angle <f> so that imposing 
that each mj is constant we get n^ 13 ' > = /A 31 ) = 0 and 
k (23) _ _ k (32) — Kq anc i therefore Eq. El yields 

/ 0 r]K 0 0 \ 

k = I —TjK 0 0 0 I . (59) 

\ 0 0 0 / 


Therefore, a rotationally invariant metamaterial around 
the Z -axis behaves as an omega-medium whose chirality 
tensor depends on a single parameter kq. This situation 
encompasses the case of a one-dimensional metamaterial 
for which e = e(Z) so that we conclude that Ea. (l59l) con¬ 
tains the most general expression for the chirality tensor 
of a one-dimensional nretamaterial, i.e. it describes one¬ 
dimensional electromagnetic chirality. 
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V. ELECTROMAGNETIC CHIRALITY OF 
ONE-DIMENSIONAL METAMATERIALS 

The description of the effective medium response, i.e. 
the evaluation of the tensors of Eas. (1551) . is achieved by 
solving Eas. flUl) for a specific dielectric distribution and 
numerical integration is usually unavoidable. However 
there is a situation where electromagnetic chirality can 
be fully discussed without resorting to numerical analysis 
and it is the case of one-dimensional metamaterials. 

Let us consider a metamaterial whose dielectric profile 
is such that e r (z) = e r (z + d), i.e. it is a periodic func¬ 
tion of the single cartesian coordinate z, whose period d 
is such that d -C A for the homogenization theory to be 
applicable. After introducing the fast spatial coordinates 
( X , Y, Z) = (x,y, z)/r] where rj = d /A is the homogeniza¬ 
tion parameter, the dielectric profile e(Z) = £ r (j]Z ) of 
Eq.® is a periodic function of period A. In Appendix 
E we show that for such dielectric distribution e(Z), the 
first of Eas. (Tl0l) can be analytically solved so that the 
effective permittivity of the first of Eas. (l3Sl) and the chi¬ 
rality tensor of Eci. (l-'16f) can be analytically evaluated and 
they are 


/ £ 
0 

0 

» ,] 


< 0 

1JK 0 

°\ 

£ 

, K = 

-r]K 0 

0 

° 

V° 

0 

e_1 ) 

\ 

K 0 

0 

°y 


(60) 

where 


averages (on the unit cell) of suitable fast-varying fields 
which, in turn, are obtained from the dielectric distribu¬ 
tion by solving a small number of electrostatic equations. 
If second order nonlocal response can be neglected, the 
medium response can be recast into the standard bian- 
isotropic form and this has allowed us to obtain a com¬ 
pact expression for the chirality tensor of the effective 
medium. We have exploited our approach to investi¬ 
gate the relation between geometrical and electromag¬ 
netic chirality by proving that the latter is shown, in ad¬ 
dition to chiral media, even by an achiral medium whose 
mirror image can be superposed onto it by means of a 
0° or 90° rotation around an axis orthogonal to the re¬ 
flection plane. We have deduced, as specific relevant ex¬ 
amples, the chirality tensor structure of two-dimensional 
and one-dimensional media and, for the formers, we have 
obtained a closed form expression for the chirality param¬ 
eter which can be evaluated from the inclusions dielectric 
profile by means of quadratures. 
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Appendix A: The electrostatic equation 


K 0 



2n 

W 




e(Z i) 

(Z 2 ) 




(61) 


The effective dielectric tensor in the first of Eqs. © co¬ 
incides with the well-known result of the zeroth-order 
homegenization theory of layered media. The structure of 
the chirality tensor in the second of Eas. (l60l) shows that 
the effective medium behaves as an omega-medium and 
this agrees with the discussion of the last paragraph, see 
Eg. (1591) . since the considered one-dimensional medium 
is in particular rotationally invariant around the Z- axis. 
It is also remarkable that the parameter kq of Ea. fliTTl) 
governing electromagnetic chirality in the present case 
can be directly evaluated from the dielectric profile by 
quadratures. 


VI. CONCLUSIONS 


Let us consider the equation 


V R • (^V R T) = $ (Al) 

whose source term $(R) has the metamaterial period¬ 
icity with vanishing spatial average $ = 0 and whose 
solution d'(R) is required to have the nretamaterial peri¬ 
odicity. Equation (lAll) can be written as id/ = >1) where 
the operator L is 


i = V R .(eV R ), (A2) 

so that solving Eq. ED amounts to the problem of in¬ 
verting the operator L on the space of functions having 
the nretamaterial periodicity. 

Note that the relation V R • (e’L SVr'Eo) = 

I 12 

e |Vr*o| + holds for any function d/o and there¬ 

fore, if 'Lo satisfies the homogeneous equation I/d'o = 0 
and it has the metamaterial periodicity we have 


To sum up we have investigated the electromagnetic 
response of a metamaterial with non-magnetic inclusions 
by a multiscale analysis and we have included spatial non¬ 
locality up to the second order. The resulting description 
is very simple and compact since the three tensors char¬ 
acterizing the effective medium response are shown to be 


[ d 3 Re|V R ’I'o| 2 = [ d? RV R • (e*SV R * 0 ) 

Jc Jc 

= I dSe%h • V r To = 0(A3) 

JdC 

where Green theorem has been used together with the 
periodicity of both e and d/g- The real and imaginary 
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parts of this equation read 


d 3 RRe(£) |V R tf 0 | = 0 


/ c 


d 3 RIm(e) |V R tf 0 f = 0 


(A4) 


and therefore, since we are considering dielectrics with¬ 
out gain for which Im(e) > 0, we obtain = 0 or 

d'o = const.. This implies that the function d'' = + H/o 

satisfies Ea. dAll) and the boundary conditions if 'L(R) 
does and therefore Eq. (ED and the periodic boundary 
conditions determines up to an arbitrary constant. 

Setting e = 1 in Eq. ED we obtain the Poisson equa¬ 
tion = $ for which the above considerations holds 

so that U/ has to be determined with the metamaterial 
periodicity up to an arbitrary constant. 


Appendix B: Tensor properties of the effective 
medium response 


A straightforward calculation shows that the relation 

= V R • [e (tf V R $ - $V R tf)] , (Bl) 

where L is the operator defined in Eo. (lA2l) . holds for any 
pair of functions W and $ so that, if both these functions 
have the metamaterial periodicity after using Green the¬ 
orem we obtain 


/ d 3 RTL$ = / d R<hLT. 

Ic Jc 


(B2) 


Let us consider a function 4'(R) which has the meta¬ 
material periodicity. The average edi'b can be rewritten 
as 

f d 3 Red l 'H = -^- [ d 3 R^e (B3) 
V Jc V Jc 

where integration by parts and the periodicity of both e 
and T have been used. Using Ea. flLfl) we obtain 


ed,JK J c d 3 T&&Lfi (B4) 

so that Eo. dB2l) with $ = fi yields 

78^ = i f d 3 KfiL . (B5) 

v Jc 

Using Ea. dBSI) with T = fj, the rank-two tensor 
Aef^^ = edifj can be written as 


a £ ; 


(eff) 


where Eq.(| 


Ae 


(eff) 

ji 


y J c d 3 RfiLfj = yf c d 3 ~RfjLfi (B6) 

]) has been used and therefore Ae < j ^ > = 


Setting T = W r j, Eq. (1B5I) becomes 

ediWrj = y J d 3 HfiLW r j 

= -y J <PB.fi[d r {efj) + Q rj ] (B7) 

where Eqs. (ED and ED have been used together with 
the property fi = 0. After integrating by parts the first 
term we get 

ediWrj = £ (<9 rfi) fj - Qrjfi■ (B8) 

Therefore the rank-three tensor = 

_ L J 1 

rje ($irfj + diW r j) becomes 

= ViQrifj - Qrjfi) (B9) 


where use has been made of the fact that fi = fj = 0. 
Setting T = v 3 Ea. dBSI) becomes 


ediVj = 


V 


d 3 RfiLvi = 


V 


d 3 RA r j£d r fi 


'c 


(BIO) 


where the second of Eas. (l33l) has been used and integra¬ 
tion by parts has been performed. Noting that Eq. ED 
and the second of Eqs. ED can be combined to yield 
sdrfi = — eS r i — p-V R A r i + Q r i we obtain 


ediVj = y J d 3 RA rj V R A r 

1 


= -eAi. 


k 2 

«, 0 


V R A ri • V R A r 


rj 


(BH) 


where the property A r j = 0 has been used and integra¬ 
tion by parts has been performed. Therefore, for the 
rank-two tensor = £ (Aij + diVj) we have 

7 = T^{JJsA r i){d s A r j). (B12) 

K o 

Setting d' = Rrsj, Ea. dBSI) becomes 


£diR rsj = y J d 3 RfiLR rs j = y J d 3 R [W sj ( ed r fi)} 

-y j d 3 R [e (. fid r W sj + 5 rs fifj)] (B13) 

where the first of Eqs. ED and Eq. ED have been used 
together with the property fi = 0 and an integration 
by parts has been performed. From Eo. d25ll we have 
sdrfi = — c5 r i — 8 r { £ fi) + Qri — LW r i which inserted in 
Eq. dB13l) , after using the property W s j = 0 and integrat¬ 
ing by parts, yields 


ediR rs j ^. 


d 3 R 


’c 


£ (d r sfifj + SriWsA-WsjLW, 


(B14) 
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With the help of this equation, the fourth-rank tensor 

fi\jlP = (fiisW r j + 5i r W S j + diR sr j + diR rs j) be- 

comes 


= ~y 

[ d 3 r( 

Jc v 


rf 

2V 


W Bi LWri + WriLWg. 


(B15) 


Using Eo. (IB2l) it is straightforward to show that P^P = 
PjisP, i.e. the tensor is symmetric with respect to the 


pair of indices ij. After symmetrizing these indices, i.e. 
setting p= \ (p^lP +Pj e isr ] ), and using Eq.® 
we obtain 

PijIP = V 2 e{-$rsfifj) 

2 

-\- — {QriW s j + QsiWrj + QrjW s i + Q s jW r i ) 

2 

-^j-£ ( fid r W s j + fidsWrj + fj8 r W s i + fj8 s Wri)- 


(B16) 


or, equivalently 

„2 


Pi/sP — -^r{QriW s j + Q s iW r j + Q r jW s i + Q s jW r i) 


V 


- ~^[fi {Prsj + Psrj) + fj (Prsi + Pat 


(B17) 


Appendix C: Effective medium electromagnetic 
chirality 

In order to show that Eas. (l45l) hold, let us first prove 
that the relation 


a 


(eff) 


1 


1 


tirk^kj — 7 CjrkK'ki 
fco ko 


(Cl) 


is satisfied. Using Eg. (PHI) we have 


(eff) 

ijr 


1/ \ (&ff) 

CC-r ' + 7T \ e irkKkj - ZjrkKki) = ® ijr 
Kq 


^irk^nmk^- 


(eff) 


^jrk^nmk^' 


ijr 

( e f f)\ 


_ e a ( e //) 6 . . 

^ t: nmq (JL nmq ^ijr 


(C2) 


which, after exploiting the relation e a b k £ c dk = dac^bd ~ 
Saddbc and the antisymmetric property = aff f \ 

becomes 


a {eff) 

ijr 


{^irk^kj tjrk^ki) — 


a (eff) 

ijr 


r (e//) 

rij 


1 


a (eff) - -e. 


v ( e //) 


jri 


2 ^nmq'~ x rimq jr 


e ijr . (C3) 


The tensor + af/f' 1 is completely 

antisymmetric since, using again the antisymmetric prop¬ 
erty it is straighforward proving that 

Fij r = — Fji r = — Fi r j = —F r ji. Since every com¬ 
pletely antisymmetric three-rank tensor in a tridimen¬ 
sional space is proportional to the Levi-Civita tensor, the 
relation Fij r = Feij r holds and the scalar F is given by 
F — Q^nmqFnmq since the relation ^-nmq^-nmq — 6 holds. 
Using the antisymmetric property we ob¬ 

tain F = SCnmq^nmq SO that 

+ a ( e fR + = p ■ =-e a {eff) e - fC4) 

^ijr ' rij ' jri 1 ^ nrn( l rimq c ijr V 

which, inserted into Ea. dC3l) . shows that Ea. dCll) is sat¬ 
isfied. 

Inserting the expression V = e, KijEj'j into the 

second of Eas. dHfl) and using the relation E' = E we 
obtain 


B' = e, ; y-K.jE’A + /i 0 H' (C5) 

which coincides with the second of Eas. flSl) . Insert¬ 
ing the expression V = ^—mjEj) into the first of 
Eas. flTfl) and using the relation E' = E we obtain 


D' = 


e i£0 


e^eo 


A e //) jq 1 / 
6 ij & j 


Jeff) 


U ^ Z-irk^kj 


(eff)-ry 1 8E'j 

s.j E j - Kki e krj dXr 


dE'j - 

8x r _ 

(C6) 


where Eq. m has been used together with the antisym¬ 
metry of the Levi-Civita tensor. Exploiting the first of 
Eas. lPHl) this equation becomes 


D' = e l£o (e[f f) E>j - ic Kkl B' k ^j (C7) 
which, using Ea. dCSI) . yields 


D' = e, 


£0 + KkiKkj ) E'j 


b^kiH k 
C 


(C8) 


which in turn coincides with the first of Eas. dTSl) . 


Appendix D: Structure of the chirality tensor 
imposed by symmetry 

In order to discuss the impact of the symmetry of the 
dielectric profile on the chirality tensor, it is essential 
to investigate preliminary the transformation rule of the 
fields fi induced by the dielectric symmetry. Differenti¬ 
ating the permittivity e(R/) with respect X[ and using 
Ea. (lt9l) we obtain 


(^R/ =Pij( d 3 e )tl- ( D1 ) 
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where Rij = RikSkj■ Let us now consider the functions 


/'(R) = /i(R') 


for which, using again Ea. (H^l) . we have 


(D2) 


[V R -(eV R /')] R 


_A_ \r(Tt'\ dMR ' y 
dX n [ [ ’ dX n 

RrnRsn [9 r (e<9 s /,;)] R ' 

[^R ‘ ( e V R /j)] R ' (D3) 


where the relation R rn R sn = S rs , stating the orthogo¬ 
nality of the improper rotation R, have been used. Com¬ 
bining Eq. (lD3l) . the first of Eas. lHUl) evaluated at R/ and 
Eq.dDD we obtain 


Vr 



-dje 


(D4) 


which structurally coincides with the first of Eas. tHOl) and 
which shows, by the uniquity of its solution, that Rijf[ = 
fj. From Eq. (ID2l) . we thus conclude that the symmetry 
of the dielectric permittivity entails for the functions /,; 
the property 

f i (R')=R ij f j (R). (D5) 


From Ea. (l46l) we note that the chirality tensor de¬ 
pends on the rank-one tensor ef m and rank-three tensor 
sfmdqfn whose geometrical properties have to be sep¬ 
arately considered. Using this equation, exploting the 
fact that the improper rotation is orthogonal (so that 
\det(R)\ = 1) and noting that the transformation in 
Eos. fl5Ul) superposes a metamaterial unit cell C onto an¬ 
other metamaterial unit cell C' we have 


efm = J 7 f d 3 R'e(R ')/ m (R') 

V Jc 

= J d 3 Re(R)R mi /i(R) = Rmisfi■ (D6) 

After differentiation, Eq. m yields 

fi ) — Rj n Rim (5 n / m ) R (D7) 

so that, reasoning as in Ea. (ID6l) we have 


efmdqfn = J d 3 R , e(R , )/ m (R / )(9 9 /„) R ' 

= y J d 3 R£(R)Rmif i (R)R qp R nj (d p f j )Yi' 

= RmiRqpRnjsfidpfj. (D8) 


Equations flU and (ID8D state that the vector ef m and 
the tensor £f m dqf n are geometrical invariants of the im¬ 
proper rotation R. 

Since R is an improper rotation (RR T = I and 
det(R) = —1), it is diagonalizable and its eigenvalues 
are Ai = e^, A 2 = e~^, A 3 = -1 (where 0 < ip < 2t r). 


In order to relate ip to the original reflection and rota¬ 
tion we note that Tr{R) = X\ + A 2 + A 3 = 2 cos ip — 1 
whereas from Eos. (|5T|) we obtain Tr(R) = RikSki = 
1 — 4 sin 2 1(rjSj) 2 so that, equating these two expres¬ 
sions, we get 


cos ip = 1 — 2 sin 2 


(D9) 


The matrix R admits three eigenvectors w^^w^ 2 ^, w®, 
such that RwA") = A„ wl"), and they satisfy the orthog¬ 
onality and completeness relations 


(n)* (m) 


= Sri 


( n )* (n) r 

W] W) = Sij. 


(DIO) 


Therefore the matrix R admits the spectral deCOmpOSi- 

Hj ~j 


tion Ra = XnW^w^* so that, after setting 


efi = VnW^, 


efidpfj = WmqnW^W^W^, 
with the help Eos. dDlOl) . Eos. dDbl) and dD8h become 


(Dll) 


Ui — A„ V n i 

tUmgp — XmX q XpWm qp 


(D12) 


where no summation is performed over the indices. 

From Eos. dD12l) it is evident that all the compo¬ 
nents of V n and W mqp vanish unless the coefficients A„ 
and X m X q X p are equal to 1. Since the allowed values 
for A n are { — 1, e 1 ^, e -1 ^} and those for X m X q X p are 
{-l,e^,e-^,-e 2i ^,-e- 2i ^,e 3 ^,-e 3 ^}, it is evident 
that the only possible cases where some of the com¬ 
ponents of V n and W mqp do not vanish are for ip = 
As a consequence these values of ip, 
corresponding to classes of original reflections and rota¬ 
tions, leads to nonvanishing chiral tensors. 

Let us now discuss the various classes of original reflec¬ 
tions and rotations corresponding to the allowed values 
of ip. In order to simplify the treatment let us consider 
the three vectors 

fl (D = fl(2) v a(3) 


il 2 ) = 


uO) = 


x ir 

u(3) x r 


\J 1 — (r • ul 3 )) 2 

— sin (|) s x r + cos (§) s 
\Jl - sin 2 (|) (f ■ s) 2 


(D13) 


which are easily seen to have real components and to 
form an orthonormal and positively oriented basis, i.e. 


. ji( m ) — A 

11 11 u nm 5 


;(«) 


x u 


( m ) - 


— 


iff) 


(D14) 
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The second of Eas. (ID14l) implies that the vectors of 
Eas. (lD13l) satisfy the relation 

CimjU$ = i e nkm . (D15) 

which is particularly useful for the developing of the 
present analysis. 

Case ip = 0. From Ea. (ID9l) we deduce that this case 
geometrically occurs both for a trivial rotation (6 = 0 ) 
and for a non trivial rotation whose unit vector is or¬ 
thogonal to the reflection unit vector (6 ^ 0 , r ■ s = 0 ). 
Since a trivial rotation is the identity transformation its 
rotation axis is arbitrary and therefore even in the first 
subcase we set r-s = 0. In both subcases, it is straightfor¬ 
ward proving that the eigenvectors of R are the vectors 
of Ecis. (ID13l) . i.e. and that, for r • s = 0, 

they reduce to 


Case ip = §,yy- From Ea. (ID9l) we deduce that 
this case geometrically occurs if the reflection and ro¬ 
tation units vectors and the rotation angle are such that 
sin 2 (|) (r • s ) 2 = 4. In this case it is straightforward 
proving that the eigenvectors are 

wW = -L(u«+m( 2 )), 

w< 2 ) = l(uW-m( 2 )), 

w (3) = u (3) . (D20) 

Since for i p = ^ we have Ai = ±i, X 2 = =R> A 3 = —1, 

from Eas. (ID12t we deduce that V\ = V 2 = V 3 = 0 and 
W mqn = 0 if it does not occur that one of the indices 
mqn is equal to 3 and the other two are both equal to 1 
or 2. Using Eas. (ID14D . (ID15I) and (ID20I) . from Eas. dDllI) 
we get after some algebra 


u (2 ^ = cos | — 1 s x r + sin 


u (3 ) = — sin s x r + cos (D16) 

Since for ip = 0 we have Ai = 1, A 2 = 1, A 3 = —1, 
from Eas. (ID12l) we deduce that V 3 = 0 and W mqn = 0 
if only one or three of the indices mqn are equal to 3. 
Using Eas. dDllI) . (ID15D and (ID16D . from Eas. dDllI) we 
get after some algebra 

< ( 3 ) (2) (2) ( 3 ) 

U\ 'Uj - U\ 'Uj 


C/mj-Z'm — V\ ^ 1 

+ V 2 , 

eimnSfmdqfn = (W 233 ~ ^322 

+ (W 112 ~ W 2 n)uf- 

+ (W331 - Wi 33 )t 4 2) t 4 3) 

+ (W 122 -W 221 )u[ 3 ) u^. (D17) 

Inserting these quantities into the definition of the chiral¬ 
ity tensor of Ea. (H51) and noting that e ?mn £/ m 3 9 /„ = 0 




efm = o, 

1 


Sfmdqfn = ~ ( W 322 ~ IU233 — W113 + W311) 


{uf ] uf> +uf ) uf ) 


+ 2 (^322 — ^233 + W113 - W311) 

{u^u^ -Mifwf) . (D21) 

Inserting these quantities into the definition of the chi¬ 
rality tensor of Eq. (l46l) we get 


Kij = T) 


» (uPu™ 


, (2) (1) N 
+ u i U j 


+K 


(a) 


(“i 


( 1 ) ( 1 ) ( 2 ) ( 2 ) 
Uj - U\ Uj 


(D22) 


where k^ (W322 - W233 - Wu 3 + W311) and 

«(“) = ^ (W322 - IU233 + W113 - W311) are two scalars 
which are mutually independent since they depends on 
four independent scalars. In matrix form the above chi¬ 
rality tensor reads 


since uW • u^ 3 ) = • u^ 3 ) = 0 we get 


K 

= v \k « (uWu^ + uP'u^) 


Kij = ?y + k^u[ 3 ^u^ 



+«(“) - u( 2 V 2 > T ) 

. (D23) 

+ «( 23 >«( 2 y. 3 > + K ( 32 U 3 y. 2)N ) 

1 J l j j 

(D18) 

Case ^ 

= From Ea.dD9l) we deduce that 


where k (13) = k 0 (V 2 +W 233 - W 332 ), 

K Oi) = k 0 {-V 2 + W 112 - W211) , k^ 

k 0 (-Vi + W331 - IU133) and k (32) = 

k 0 (Ui + W 122 — W 221 ) are four scalars which are 
mutually independent since they depends on ten in¬ 
dependent scalars. In matrix form the above chirality 
tensor reads 


n = 77 


(> 3 )u«u 


( 3 ) T + /«(31)u(3)u(l) T 


this case geometrically occurs if the reflection and rota¬ 
tion units vectors and the rotation angle are such that 
sin 2 (|) (r • s ) 2 = |. Since for ip = ^^ we have 
Ai = \ (- -1 ± iy/ 3), A 2 = 5 (—1 T iy/3 ), A 3 = —1, from 
Eas. dD12l) we deduce that V\ = V 2 = V 3 = 0 and 
W m qn = 0 if it does not occur that m = q = n = 1 
orm = q = n = 2 (i.e. only the components Wm and 
IU 222 survive). From Eas. dDllI) we directly get 


+ k^u^u^ t + k^u^u^) . (D19) 


c 27717 


e/m = 0 , 


lefmdqfn = 0 


(D24) 
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so that from the definition of the chirality tensor of 
Eq.® we obtain 


Kij = 0, (D25) 

i.e. the chirality tensor vanishes. 


Appendix E: Effective medium response of 
one-dimensional dielectric media 


Let us consider a metamaterial whose dielectric profile 
is e(Z), i.e. it is dependent on a single cartesian coordi¬ 
nate. The dielectric profile e(Z) and its inverse [e(Z )] _1 
can be represented through the Fourier series 

+ OO 

e(Z) = e ' inkoZa n, 

n=—oo 

+oo 

\e{Z)]- 1 = Y e inkoZ bn (El) 


where ko = 2 tt/X and 

a n = \ [ X dZe~ ink ° z e(Z), 

A Jo 

HZ)\ 


bn — 


dZe 


(E2) 


are their Fourier coefficients. Due to the generality of 
the Fourier series representation, the situation we are 
here considering allows us to investigate even discontinu¬ 
ous dielectric distributions, a relevant situation occurring 
when the metamaterial is a stratified layered structure, 
without explicitly resorting to matching condition on the 
surfaces of the inclusions. Note also that 


a 0 = e, 

b 0 = F 1 . (E3) 


For the one-dimensional dielectric distribution e(Z), 
the functions fj satisfying the first of Eqs. ® are evi¬ 
dently dependent solely on Z, i.e. fj = fj{Z), so that 
such equations yields 


d_ 

dZ 

_d_ 

dZ 



0, j = 1,2, 
de 

~dZ' 


(E4) 


which have to be solved with the requirements that the 
functions fj have the same periodicity as s(Z) and have 
vanishing average fj =0. After an integration Ecis. (IE4l) 
become 


dfi 

dZ 

dfz 

dZ 


Cl 

£ ’ 

-1 


3 = 1,2, 


Cz 

£ 


(E5) 


where Ci, C 2 , C 3 are constants. Averaging Eas. (IE5l) over 
the metamaterial period, using Eas. (lE3|) and noting that 
the average of their left hand sides are zero since they are 
derivatives of periodic functions we get 


0 = Cjb 0 , j = 1,2, 

0 = -I + C 360 . (E 6 ) 

which, assuming bo 0, imply that C\ = C 2 = 0 and 
C 3 = I/ 60 ■ Accordingly, from the first of Eas. (IE5l) we 
obtain that /1 and fi are constants so that, requiring 
their average to be zero, we get 

fi = h = 0. (E7) 

Using the Fourier espansion of e _1 of the second of 
Eqs. (EH), the second of Eas. dESI) becomes 

§- £ < E8 > 

n=— oo,n^0 

which, after an integration and the requirement for /3 to 
have vanishing average, yields 


h 


+ OO 

E 


ginfeo Z 


n——oo,n^0 


inkobo 


(E9) 


Combining Eas. dE7l) and dE9l) we obtain for the functions 

fj 


fj=*. 


13 


+ OO 

'y ^ ginkoZ 


n=—oo,n^0 


inkobo 


(E10) 


which are fully determined by the dielectric profile e(Z) 
and they are unique as a consequence of the chosen con¬ 
straints (periodicity and vanishing average). 

Using Ea. dElOl) . the first of Eas. dMl) becomes 


( +00 \ / +00 , 

Y eimkoZ am) [Sij+SxSjs Y einkoZ Y 

m— — 00 / \ n=—oo,n^0 ' 

(Ell) 

which, inserted into the first of Eas. (l38l) and neglecting 
the rj 2 contribution, yields 


^ +00 

= §. jCl0 + 8 i 3 Sj 3 — Y a-nbn■ (E12) 

n——oo,n^0 


Averaging the relation ee 1 = 1 and using Eqs. m 
the relation Jf,n=-oo a -nb n = 1 easily follows so that 
J2n=-oo,n?0 a -nbn = 1 - a 0 b 0 and Eq. dE12|) becomes 


Je//) 

ij 


{SnS. 


'll 


Si25j2)ao + 5i3Sj3 — . 

00 


or in matrix form 


£ (e//) 


/ e 0 
0 £ 

^0 0 



(E13) 


(E14) 
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where Eas. (IE3l) have been used. Inserting Eas. (IE2l) into Ea. (IE17|) we obtain 

In order to evaluate the chirality tensor we note that 
Eos. (IE 101) and the first of Eas. dEll) yield 


+ OO 


£fm — fim3 ^ ^ 


s=—oo,s ^0 


CL —§ b§ 

ikobos ’ 


efmdqfn = S m3 6 q3 S n3 ef 3 ^ 

so that, from Eq. fl^l) we obtain 

Kij — C-ij^TjHLQ 

where 


+ °° 7 

i v > d— s b s 

Ko = T - 

bo ^ , s 

s=— 00,57=1) 


In matrix form, Eo. (IE16|) reads 

0 rjno 0 

k = | —T]t Co 0 0 

0 0 0 


(E15) 


(E16) 


(E17) 


K 0 = 


boX 2 


rX dZl r dZ ^ 


’OA J o 


/o 


-s( Z2 ) 


+oo 

E : 

s=—oo,s ^0 


pi27rs (£i^2) 


(E19) 

After noting that the even contribution in the numerator 
within the series does not contribute and using the rela- 
tion XEn sln( f si) = 7T [-£ + isign(0] valid for |£| < 1, 
Ea. (IE20l) yields 


k 0 = 


‘I /V' 1 - E(Zi 


2 S1 § n 


dZ,2 ' F (Z 1 
0 Jo £ \ Z/ 2) 

Z \ — Z 2 


Z \ — Z 1 


A 


(E20) 


(E18) where the first of Eas. (IE3l) has been used. 
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